Upward muon signals at neutrino detectors as a probe of dark matter properties 
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We study upward muon flux at neutrino detectors sucli as Super-Kamiokande resulting from 
high-energy neutrinos produced by the dark matter annihilation/decay at the Galactic center. In 
particular, we distinguish showering and non-showering muons as their energy loss processes inside 
the detector, and show that this information is useful for discriminating dark matter models. 
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I. INTRODUCTION 

Recent observations of cosmic ray positron/electron 
fluxes by the PAMELA satehite [l] and ATIC bal- 
loon experiments [2] have shown excess of events com- 
pared with expected astrophysical backgrounds. PPB- 
BETS [3], which is another balloon experiment, also re- 
ported cosmic ray electron flux consistent with ATIC 
results. Although more recent results of Fermi satel- 
lite and HESS telescope [5j] on the electron flux did 
not conflrm such a sharp excess observed by ATIC, the 
positron excess observed by PAMELA still remains vi- 
able. While these anomalous positron/electron fluxes 
may be explained by some nearby high-energy astrophys- 
ical sources [6, 7], they can also be interpreted as evidence 
of dark matter annihilation/decay [8]. If the recent Fermi 
data is fltted by dark matter annihilation/decay, a heavy 
dark matter with mass of order a few TeV is favored [9J. 
In the case of annihilation, the required annihilation cross 
section for reproducing the PAMELA anomaly is around 
{(jv) ~ 10~^^ cm^s~^, which is orders of magnitude larger 
than the standard value, 3 x 10~^^ cm^s~^. In order to 
flll this gap, one needs to invoke nonthermal dark matter 
production scenario [TO, El, Ell or Sommerfeld enhance- 
ment mechanism [l3| giving velocity-dependent annihila- 
tion rate. In the case of decay, the required decay rate is 
about r - 10-2^ s-^ 

In order to conflrm that the positron/electron excesses 
are caused by dark matter annihilation/decay, cross 
checks by other observations are necessary. Gamma- 
rays anti-protons [l5[, anti-deuterons [l^ and ra- 
dio emission associated with the same annihilation/decay 
processes [17] may be used for this purpose. At the 
present stage, however, no clear excess was found on 
those observations, and they give constraints on anni- 
hilating/decaying dark matter models. Big-bang nucle- 
osynthesis [18|, ll9| and the reionization |20] also give 
stringent constraints on annihilating dark matter mod- 
els. 

It was also pointed out that dark matter anni- 
hilation/decay which explains the PAMELA anomaly 
may predict potentially large neutrino flux as well [21, 
[22I . which are constrained from observations of Super- 
Kamiokande (SK) toward the direction of the Galac- 



tic center [23|, |2^. Dark matter annihilation (decay) 
yields high-energy neutrinos at the Galactic center where 
the number density is enhanced, and they reach to the 
Earth without energy loss. During the propagation, the 
neutrino oscillation effects mix three flavors (i^g, ^/i, ^r)- 
Among them, some fraction of muon neutrinos are con- 
verted into high-energy muons due to the interaction with 
matter inside the Earth, and if the interaction point is 
within around 1 km from the detector, produced muons 
penetrate the detector leaving characteristic signatures 
like the Cherenkov light from a muon itself and/or from 
electromagnetic shower produced by it [25|, [lO] • Since the 
direction of a primary neutrino is conserved even after it 
is converted into a muon, SK limit on the muon flux from 
the Galactic center direction provides useful constraint 
on dark matter annihilation/decay models [0, 

In this paper we further study the neutrino flux from 
the Galactic center generated by dark matter annihila- 
tion/decay and how they are detected at the neutrino de- 
tector. Since both conversion rate of neutrino into muon 
and muon energy loss processes in the Earth matter de- 
pend on its energy, the signal of the high-energy neutrino 
has energy dependence. In particular, we distinguish 
showering and non- showering muon events at SK, follow- 
ing the criterion given in Ref. [2J]. As is expected, higher 
energy muons are likely to lose energy by radiative losses, 
developing electromagnetic shower inside the detector. 
Interestingly, the probability for causing shower events 
becomes signiflcant around the muon energy ~ 1 TeV, 
which is a natural scale of dark matter mass. Therefore 
such a separation, showering/non-showering events, has 
rich information on the initial neutrino energy spectrum 
and we show that it is a possible method to discriminate 
dark matter properties.^ This analysis can be applied to 
future neutrino detectors, such as KMSNeT, for search- 
ing for dark matter, or even the present on-going exper- 
iments, such as SK, may be able to improve constraints 
on dark matter properties. 

In Sec. iniwe provide formulae to evaluate the upward 
showering and non-showering muon fluxes from the pri- 



^ The importance of a separation of showering and non-showering 
events was also pointed out in Ref. [9|]. 
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mary neutrino flux generated by dark matter annihila- 
tion/decay at the Galactic center. Applications to some 
dark matter annihilation/decay models is performed in 
Sec. mil Sec. IIVI is devoted to conclusions. 



II. MUON FLUX FROM 
ANNIHILATING/DECAYING DARK MATTER 

In this section we develop a formalism to evaluate 
the upward muon flux starting from the initial neu- 
trino energy spectrum produced by dark matter anni- 
hilation/decay. 



A. Primary neutrino flux 

We assume that DM particle, denoted by with mass 
m annihilates (decays) into some flnal state F with an- 
nihilation cross section (decay rate) {av)^ (I^f)- The 
primary neutrino flux from the Galactic center is given 
by 



dN 



for the case of annihilation, and 
for the case of decay, with 




(J2)QAf^, (1) 



(Ji)nAn, (2) 



dN 



dN 



dE 



(3) 



E=E^ 



where i = e,/i,r, and P^iu^ denotes the probability that 
the Ui at the production is observed as V/j, at the Earth 
due to the effect of neutrino oscillation and given as 
P^^^^ = 0.22, P^^^^ = 0.39, and P^^^^ = 0.39. Here, 
Rq = 8.5 kpc and = 0.3 GeVcm"^ are the distance 
of the solar system from the Galactic center and the local 
dark matter density near the solar system, respectively. 
F collectively denotes the primary annihilation/decay 

modes (e.g., etc.), and dNp^'^/dE represents the 

Ui spectrum arising from the flnal state F which is sim- 
ulated by PYTHIA package [2^. The dependence on 
the dark matter halo density proflle is contained in the 
remaining factor {Jn)n^ deflned by 



{Jn)Q 



r dQ r 
J ^L. 



dm (p{i)y 

Rq \Pq J 



(4) 



where /(?/;) is the distance from us along the direction 
which is the cone-half angle from the Galactic cen- 
ter within the range < t/j < V^max, and AQ{= 27r(l — 
cos^/^max)) IS the solid angle over which the neutrino flux 
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TABLE I: Typical values of J-factor for NFW and isothermal 
profiles for ^max = 5°, 10°, 15°, 20° and 25°. 



is averaged. Dark matter halo proflle is characterized as 
follows: 



p{r) 



Po 



(r/i?)7[l + (r/i?)«](/3-7)/«' 



(5) 



where r is the distance from Galactic center. Here 
we consider two typical density proflles: Navarro- Frenk- 
White (NFW) proflle [29] (a = 1, /3 = 3, 7 = 1, = 
20 kpc) and cored isothermal proflle (a = 2,/? = 2,7 = 
0, = 2.8 kpc). Typical values of {Jn)nAQ for these 
proflles are summarized in Table [H As is expected, dark 
matter halo with large core |30| predicts smaller neu- 
trino flux, but the dependence on the halo proflle is not 
so large compared with gamma-ray flux from the Galac- 
tic center [ij , as long as the angular resolution is larger 
than ~ 5°. Thus following discussion does not much de- 
pend on the halo proflle and it is possible that neutrino 
signals are detectable while the gamma-ray constraints 
are satisfled. 



B. Upward muon flux 

High-energy neutrinos produced by the dark matter 
annihilation/decay at the Galactic center come in the 
Earth and some of them interact with matter inside the 
Earth, generating high-energy muons. Since muons have 
electric charges, they lose their energy during the prop- 
agation inside the Earth. Typical propagation distance 
of high-energy muons with energy ~ 1 TeV is around 
1 km, as will be shown. Therefore, muons appearing only 
within 1 km around the detector can reach the detector, 
leaving their characteristic signals. The search strategy is 
to flnd upward muons in order to avoid cosmic-ray muon 
background. 

Taking these processes into account, showering muon 
flux from the Galactic center within cone-half angle is 
given by 



(6) 



where dFj,^/dEj,^ is the parent differential neutrino flux 
evaluated in the previous subsection, Ej^^ is the parent 
neutrino energy, f{Ej^^) denotes the probability that a 
neutrino with energy Ej^^ is converted into a muon with 
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FIG. 1: Showering fraction e(Ejy ) as a function of initial 
neutrino energy Eu^ given in Ref. [23 ]. 



observable energy range, and e{Ei,^) is fraction of the 
showering muon in the total observed muon events with 
parent neutrino energy Ej^^ . 

The last ingredient, e{Ey^)^ is read off from Fig. 5 of 
Ref. [23]. By taking the ratio between black and black 
plus red lines in that figure, we obtain a following fitting 
formula. 



bx - 



(7) 



where x = \og^Q{Ej,JlGeV) and a = -0.0491232,6 = 
0.059575, c = -0.0162841, and d = 0.0103345. (Note 
that this differs from the quantity e{E^) used in 
Ref. The result is shown in Fig. [H Note that 

this formula is valid for 10 GeV < E^^ < 10 TeV. If we 
want to know total muon events, including non-showering 
events, we have only to replace e{Ej^^) with 1: 



^(total) 



dE, 



dF,, 



dE, 



(8) 



In order to evaluate f{Ej,^)^ the muon-nucleon scat- 
tering cross section and the muon energy loss in the rock 
surrounding the SK must be known. It is given by 



f(T^ \ Z^"'' ^ZT^ ( ^^^^ {v^.)p^^^{^l)X (rock) 



(9) 



muon with energy given by [31^^ 



dE, 



^ ' 0.29 + 0.21-1- 1 . (11) 



E?. 



The mass density of the standard rock is 2.65 gem ^ 
[32[, and hence rip^^^-' = nn°^^^ = 2.65A^a/2 cm~^ where 
Na = 6.022 X 10^^ is the Avogadro's number. R{E^, Eth) 
is the distance which muons with energy E^ can travel 
until their energies decrease to the threshold energy £^th 
in the rock. It is obtained by solving the energy loss 
equation. 



dX 



(12) 



where X represents the distance in units of gcm~^. The 
first term a{E^) is from ionization loss, and it only 
logarithmically depends on E^. The second term de- 
notes other radiative energy loss effects due to bremm- 
strahlung and pair creation, which may cause electromag- 
netic shower inside the detector. The precise formulae for 
those energy losses are given in Ref. [32]. In Fig. El we 
show the propagation distance of muons in the rock. For 
reference, we also have shown the result when only the 
ionization loss is taken into account. For convenience, we 
provide fitting formula for R{E^^Eth = 10 GeV) in the 
rock, which is valid for 10 GeV < < 10 TeV, as 



R{E^,EtY, = 10 GeV) = 10 



a-\-hy-\-cy 



[kr 



(13) 



where y = \og^Q{E^/lGey) and a = -3.29186,6 = 
1.52594, and c = -0.147224. 

Notice that both the interaction cross section of neu- 
trinos with nucleon and propagation distance of muons 
are roughly proportional to Ej^^{^ E^). Therefore heav- 
ier dark matter has a tendency to provide more muon 
signals. This situation is quite different from other indi- 
rect signals, such as gamma-rays or positrons. For fixed 
annihilation cross section or decay rate, smaller fluxes 
of gamma-ray or positron will be predicted for heavier 
dark matter, since the number density of dark matter is 
inversely proportional to the dark matter mass. In the 
case of neutrino-induced muon signals, however, higher 
energy neutrinos more likely produce muon signals, so 
that it compensates the suppression from smaller num- 
ber density of dark matter for heavier mass [21i] . 



where n^p^^^^ is the number density of proton (neutron) in 



the rock and da, 



I dE^ is the scattering cross 



section of a neutrino with proton (neutron) to create a 



^ This includes both processes such as v^'p fiX and u^p ftX. 
Thus primary neutrino flux in Eq. (|6]) should only count neutri- 
nos, not anti-neutrinos (assuming that primary flux of neutrinos 
and anti-neutrinos are equal). 
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FIG. 2: Propagation distance of muon in the standard rock 
as a function of muon energy, when all energy loss processes 
are taken into account (solid line) and only ionization loss is 
considered (dashed line). 



C. Atmospheric background 

Before going into predictions of dark matter models, 
we evaluate the background muon events coming from 
atmospheric neutrinos. Since we are interested in the 
upward muon events, we do not need to worry about 
background from atmospheric muons. Detailed calcu- 
lations for the atmospheric neutrino spectrum are per- 



formed in Refs. [3j 
flux given in Ref. 




events are evaluated as 

; 0.34 X IQ- 



Using the atmospheric neutrino 
the background upward muon 
2.8x10-1^ cm-^s-^ 

~^s~^ from the cone 



-15 



(corresponding to N^^^^l'^^"^ 



sr 



and iv(*°-<''-'«G) ^ ^ 4 



and iv^t"!'*''''^^^ 
half angle V^max 
1.2 X IQ-i^ cm 

10~^^ cm~^s~^sr~^, respectively). If dark matter predic- 
tion is larger than or comparable to those values, dark 
matter detection by neutrino-induced muon signals is 
promising. 



III. DISCRIMINATING DARK MATTER 
MODELS 

Now let us evaluate the upward showering/non- 
showering muon flux following the procedure given in the 
previous section. We consider dark matter mainly anni- 
hilating (decaying) into r+r~, ZZ, and 
^i^i = e,/i,r). Following results are not changed so 
much between and r+r", and hence we call these 
models as "leptonic annihilation (decay)". Similarly, an- 
nihilating (decaying) into and ZZ lead to almost 
the same results, and we call those cases as "hadronic 
annihilation (decay)". Taking account of neutrino oscil- 
lation effects, annihilating (decaying) into z^^z^^ and Ur^r 



yield completely the same muon flux, but the case of z/e^e 
predicts slightly smaller flux. In the following we call the 
former case as "neutrino annihilation (decay)", since the 
characteristic behaviors are essentially the same. 

Resulting total and showering muon fluxes are shown 
in Figs. [3]|5l for the above mentioned three cases: neu- 
trino, leptonic, and hadronic annihilation (decay). In 
these flgures, the vertical axis is {av){J2)nAQ for the 
case of annihilation (top panels), and T{Ji)qAQ for the 
case of decay (bottom panels). Contours of Nj^^^^P = 

(0.1, 1, 10) X 10-15 cm-2s-i and AT^?;""^ = (0.1, 1, 10) x 
IQ-i^ cm-^s-i are shown by solid and dashed lines, re- 
spectively. The shaded region is excluded at 90% C.L. 
by the observation of SK toward the direction of Galac- 
tic center with cone half angle 5° [23]. If other values 
of cone angle are chosen, the constrained region becomes 
widened or narrowed. 

As is easily seen from these flgures, in the case of an- 
nihilation, the total muon flux remains roughly constant 
even if dark matter mass is increased for flxed annihi- 
lation cross section. This is because the annihilation 
rate is proportional to m"^ reflecting the fact that the 
rate is proportional to the square of dark matter num- 
ber density, but the muon production cross section inside 
the Earth and the muon propagation distance are both 
roughly proportional to m, which cancels the m depen- 
dence. However, the muon propagation distance deviates 
from the simple scaling R{E^) ex when the radiative 
energy loss becomes important at Ej^ > 500 GeV (see 
Fig. [2]). Thus for heavy dark matter with m > 500 GeV, 
the predicted total muon flux shows mass dependence. 
On the other hand, the showering muon is caused by 
high-energy component of the primary neutrino flux, and 
hence its mass dependence is quite different from the 
total muon flux, although the flux itself is suppressed. 
However, the showering muon flux prefers heavier dark 
matter mass, and atmospheric background flux is signif- 
icantly reduced for higher energy region. Thus we may 
have a chance to detect muon signals and observe sepa- 
rately non-showering and showering muons. 

Mass scaling in the case of decaying dark matter can 
also be recognized in a similar way. Since the decay rate 
per unit volume is only proportional to the dark matter 
number density, the total muon flux is proportional to m 
for flxed decay rate. This feature can be read off from 
the flgures. 

Fig. [6] shows the ratio of the shower/total muon flux 
for various flnal states: z^z^, r+r~, and 
ZZ. This ratio is not sensitive to the neutrino flavors: 
and v^Vr yield a common ratio. It is seen 
that the shower/total muon ratio is quite sensitive to the 
dark matter mass irrespective to the dark matter annihi- 
lation/decay modes. Since this ratio expected from the 
atmospheric neutrino background is about 0.12, the dark 
matter signals predict clearly large shower/total muon ra- 
tio. Thus searching for muons signals with distinguishing 
showering and non-showering muons from particular di- 
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FIG. 3: Contours of n'^T"'- 



and N 



(0.1,1, 

(shower) ^ (0.1,1,^10) X 10-^^ Cr 



10) X 10"^^ cm"^s"^ 

i^^s""*^ are shown by 
solid and dashed hnes, respectively. The horizontal axis is 
dark matter mass and the vertical axis is {av){J2)nAQ for 
the case of annihilation (top), and r(Ji)nAQ for the case of 
decay (bottom). Here dark matter is assumed to annihilate 
(decay) into z^^z>/^. Results in the case of ly^i^T are the same. In 
the case of i^e^^e, both fluxes are slightly reduced by an amount 
of factor 2. The shaded region is excluded at 90% C.L. by the 
observation of SK toward the direction of Galactic center with 
cone half angle 5°. 



rection may be useful for discriminating the signal from 
the background and constraining/extracting properties 
of dark matter. 



IV. CONCLUSIONS 

In this paper, we have studied how dark- matter origi- 
nated neutrinos can be detected and how can we extract 
useful information on the dark matter properties, mo- 
tivated by recent observations of anomalous cosmic-ray 
positron/electron fluxes. 

For this purpose, we have focused on the separation be- 
tween showering and non-showering upward muon fluxes 
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FIG. 4: Same as Fig. O but for annihilating (decaying) into 
In the case of t^t~, muon fluxes are slightly smaller 
by an amount of ~ 10%. 



inside the neutrino detector. Since muons with energy 
larger than about 500 GeV lose their energy dominantly 
by radiative processes, such as bremmstrahlung and pair 
creation processes, searching for electromagnetic shower 
is a useful way to discriminate signals of high-energy neu- 
trinos from atmospheric neutrino background. Moreover, 
since the energy dependence of the total muon flux and 
showering muon flux are different from each other, mea- 
suring the ratio between them gives us information on 
the dark matter mass. A good point of muon search 
compared with gamma-ray and positron/electron flux is 
that the muon flux is enhanced for heavier dark mat- 
ter due to the enhancement of the muon production and 
propagation distance inside the Earth. Thus our method 
has a particular benefit for rather heavy dark matter 
case, as is favored from the recent observations of the 
cosmic ray electron flux by the Fermi satellite. The pre- 
vious bound on the neutrino flux from the Galactic center 
given in Ref. [23] comes from data of SK-I. Since the Su- 
perKamiokande accumulates more data now, the bound 
is expected to be improved soon. 

Finally we mention possible applications to the on- 
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FIG. 5: Same as Fig. O but for annihilating (decaying) into 
VK^VK". Results in the case of ZZ are almost the same. 



as dwarf galaxies may be used as a source of high-energy 

neutrinos. Future megaton scale water tank detector, 

Hyper-Kamiokande, and kilometer size detector, such as 

KMSNeT, is expected to significantly improve current 

sensitivities of SK and may be ideal experiments to search 

for the dark matter signatures by looking at the Galactic 

center. Detailed estimates on the sensitivities in these 

experiments are beyond the scope of this paper, and we 

leave this issue for future work. 
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FIG. 6: Ratio of the showering muon to the total muon fluxes 
as a function of dark matter mass in the dark matter an- 
nihilation. In the case of decay, the horizontal axis corre- 
sponds to half of the dark matter mass. The label shows 
the dominant mode of the dark matter annihilation/decay: 
zyzy,/i+/i-,r+r", V^+H^", and ZZ. 



going or future neutrino detectors. The IceCube detector 
cannot see the upward muon signals from the Galactic 
center, since it is located at the South pole. A planned 
extension, called IceCube DeepCore [35|, may have a po- 
tential to search Galactic center by looking at downward 
muons if atmospheric muons can be removed to the level 
of expected signals. Although the sensitivity of the Deep- 
Core on the muon flux from the Galactic center is diffi- 
cult to estimate at the present stage, other targets such 
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